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Abstract 

Dijet production in deep inelastic ep scattering is investigated in the region of low 
values of the Bjorken-variable x (10~ 4 < x < 10~ 2 ) and low photon virtualities Q 2 (5 < 
Q 2 < 100 GeV 2 ). The measured dijet cross sections are compared with perturbative QCD 
calculations in next-to-leading order. For most dijet variables studied, these calculations 
can provide a reasonable description of the data over the full phase space region covered, 
including the region of very low x. However, large discrepancies are observed for events 
with small separation in azimuth between the two highest transverse momentum jets. This 
region of phase space is described better by predictions based on the CCFM evolution 
equation, which incorporates k t factorized unintegrated partem distributions. A reasonable 
description is also obtained using the Color Dipole Model or models incorporating virtual 
photon structure. 
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1 Introduction 



Dijet production in deep inelastic lepton-proton scattering (DIS) provides an important test- 
ing ground for Quantum Chromodynamics (QCD). At HERA, data are collected over a large 
range of the negative four-momentum transfer squared, Q 2 , the Bjorken-variable x and the trans- 
verse energy, E T , of the observed jets. HERA dijet data may be used to gain insight into the 
dynamics of the parton cascade exchanged in Iow-jc lepton-proton interactions. Since in this 
region of phase space photon-gluon fusion (Figured) is the dominant underlying process for 
dijet production, such measurements open the possibility of studying the unintegrated gluon 
distribution, first introduced in JTJE1- 

In leading order (LO), i.e. 0(a s ), dijet production in DIS is described by the boson-gluon 
fusion and QCD-Compton processes (Figure[T^, b). The cross section depends on the fractional 
momentum x of the incoming parton, where the probability distribution of x is given by the 
parton density functions (PDFs) of the proton. The evolution of the PDFs with the factorization 
scale, fil, is generally described by the DGLAP equations [3 1. To leading logarithmic accuracy, 
this is equivalent to the exchange of a parton cascade, with the exchanged partons strongly 
ordered in virtuality up to Q 2 . For low x this becomes approximately an ordering in k t , the 
transverse momentum of the partons in the cascade (Fig. [TJ;). This paradigm has been highly 
successful in the description of jet production at HERA at large values of Q 2 or E 2 T of the 
jet 00012. 

The DGLAP approximation is expected to breakdown at low x, as it only resums leading 
logarithms in Q 2 and neglects contributions from log 1 fx terms, which are present in the full 
perturbative expansion. This breakdown may have been observed in forward jet and forward 
particle production at HERA lEM fTWTl . 

Several theoretical approaches exist which account for \ow-x effects not incorporated into 
the standard DGLAP approach. At very low values of x it is believed that the theoretically 
most appropriate description is given by the BFKL evolution equations which resum 

large logarithms of l/x up to all orders. The BFKL resummation imposes no restriction on the 
ordering of the transverse momenta within the parton cascade. Thus off-shell matrix elements 
have to be used together with an unintegrated gluon distribution function, f(x,p. 2 ,k t ), which 
depends on the gluon transverse momentum k, as well as x and a hard scale fif. A promising 
approach to parton evolution at low and larger values of x is given by the CCFM |[P3l evolution 
equation, which, by means of angular-ordered parton emission, is equivalent to the BFKL ansatz 
for x — » 0, while reproducing the DGLAP equations at large x. 

Experimentally, deviations from the DGLAP approach may best be observed by selecting 
events in a phase space region where the main assumption, the strong ordering in k t of the ex- 
changed partons in the cascade, is no longer expected to be fulfilled. This is the case at low x. 
Parton emission along the exchanged gluon ladder (Fig. [TJ:) increases with decreasing x. This 
may lead to large transverse momenta of the partons entering the hard scattering process, such 
that in the hadronic (photon-proton) center-of-mass system (cms) the two partons produced in 
the hard scattering process (Fig[T^,b) are no longer balanced in transverse momentum. For the 
final state studied here, one then expects an excess of events in which the two hardest jets are no 
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longer back-to-back in azimuth. Such configurations are not included in DGLAP-based calcu- 
lations, necessitating the inclusion of additional contributions when calculating cross sections 
at low x. 

An alternative approach to modelling additional contributions ifbH due to non-fc r ordered 
parton cascades is given by the concept of virtual photon structure. This approach mimics 
higher order QCD effects at low x by introducing a second fc r ordered parton cascade on the 
resolved photon side, evolving according to the DGLAP formalism. This resolved contribution 
is expected to contribute for squared transverse jet energies, Ej., greater than Q 2 , which is the 
case for most of the phase space of the present analysis. Virtual photon structure is expected to 
be suppressed with increasing Q 2 . Leading-order QCD models which include the effects of a 
resolved component to the virtual photon have been successful in describing dijet production at 

low Q 2 ca. 

The aim of this paper is to provide new data in order to identify those regions of phase 
space in which next-to-leading order (NLO) DGLAP-based QCD calculations are able to cor- 
rectly describe the underlying dynamics of the exchanged parton cascade and those in which 
the measurements deviate from the DGLAP based predictions. Where deviations are observed, 
comparisons of the data with other QCD models are performed. Dijet production at low x and 
low Q 2 is an appropriate tool for this purpose as the jet topology reflects the dynamics of the 
parton cascade lfT6lfT71fTH1 . Therefore, dijet cross sections are measured multi-differentially 
as a function of observables particularly sensitive to Iow-jc dynamics, considerably extending 
an earlier analysis |[T9l in terms of the observables studied, the kinematic reach and statistical 
precision. 



2 Experimental Environment 

The measurement presented is based on data collected with the HI detector at HERA during 
the years 1996 and 1997. During this period the HERA collider was operated with positrons 1 of 
27.6 GeV energy and protons with energy of 820 GeV The data set used corresponds to a total 
integrated luminosity of 21 pb _1 . 

The HI detector consists of a number of sub-detectors [ 20 1 providing complementary and 
redundant measurements of various aspects of the final state of high energy electron-proton 
collisions. The detector components which are most important for this analysis are the backward 
calorimeter, SpaCal [21 1, together with the backward drift chamber, BDC ll2"2"l . for identifying 
the scattered electron, and the Liquid- Argon (LAr) calorimeter [ 23 ] for the measurement of 
the hadronic final state. The central tracking system is used for the determination of the event 
vertex and to improve the hadronic energy measurement by the LAr calorimeter [24]. 

The SpaCal is a lead/scintillating-fiber calorimeter covering polar angles 2 in the range 1 53° < 
< 177.5°. Its electromagnetic part has a depth of 28 radiation lengths and provides an energy 

'in this paper we refer to the incident and scattered lepton as "electron". 

2 The z axis of the right-handed coordinate system used by HI is defined to lie along the direction of the proton 
beam with the origin at the nominal ep interaction vertex. 
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resolution of <t e /E « 0.07/ V£[GeV]©0.01 [25 1. Remaining leakage of electromagnetic show- 
ers and energy depositions of hadrons are measured in the hadronic part of the SpaCal. The 
accuracy of the polar angle measurement of the scattered electron, using the vertex position and 
the BDC (156° < 9 < 175°) is 0.5 mrad ESI . 

The LAr calorimeter covers the angular region 4° < 6 < 154°. Its total depth varies between 
4.5 and 8 interaction lengths, depending on the polar angle. It has an energy resolution of 
cr E /E « 0.50/ V£[GeV] © 0.02 for charged pions E71l . The LAr calorimeter surrounds the 
central tracking system, which consists of multi-wire proportional chambers and drift chambers, 
providing measurements of charged particles with polar angles of 15° < 9 < 165°. 

The SpaCal and LAr calorimeters are surrounded by a superconducting solenoid, which 
provides a uniform field of 1.15 T parallel to the beam axis in the region of the tracking system, 
allowing track momentum measurements. 

The luminosity is determined from the rate of the Bethe-Heitler process {ep — » epy). The 
luminosity monitor consists of an electron tagger and a photon detector, both located down- 
stream of the interaction point in the electron beam direction. 

3 Selection Criteria 

The analysis is based on a sample of DIS events with a clear multi-jet topology of the hadronic 
final state. The events are characterized by an electron scattered into the backward calorimeter, 
SpaCal, and at least two jets within the acceptance of the LAr calorimeter. They are triggered 
by demanding a localized energy deposition in the SpaCal and by track requirements, which 
result in a trigger efficiency of (97.3 + 0.1)% [T2"8l. 

The scattered electron is identified as the cluster of highest energy, E e > 9 GeV, in the 
electromagnetic part of the SpaCal. In order to select well identified electromagnetic showers, 
a cut of 3.5 cm is applied on the energy weighted radius of the selected cluster [28 1. The energy 
in the hadronic part of the SpaCal within a radius of 15 cm of the shower axis is required to be 
less than 0.5 GeV. Moreover, an electron candidate must be associated with a track segment in 
the BDC. 

The inclusive event kinematics are derived from the energy and polar angle measurements 
of the electron candidate. The kinematic range of the analysis is restricted to the low-Q 2 , low-* 
region, 5 < Q 2 < 100 GeV 2 and 10" 4 < x < 10~ 2 . In addition, the inelasticity y = Q 2 /xs 
is restricted to 0.1 < y < 0.7, where V* ls me ep center-of-mass energy. In the acceptance 
region of the SpaCal, the restriction y < 0.7 always corresponds to the requirement E e > 
9 GeV on the energy of the scattered electron. The requirement y > 0.1 ensures a large central 
track multiplicity and hence the accurate reconstruction of the event vertex, which is required 
to lie within |z vtx | < 35 cm. The restriction of the z/-range also reduces the effects of QED 
bremsstrahlung. 

The requirement of E e > 9 GeV suppresses background from photoproduction processes in 
which the scattered electron escapes through the beam-pipe but an electron signal is mimicked 
by a particle from the hadronic final state. This background is further reduced by demanding 
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35 < Yii(Ei - p z j) < 70 GeV. Here the sum runs over the energies and momenta of all final state 
particles including the scattered electron. For fully reconstructed events, energy and momentum 
conservation implies that Yii(Ei- p z j) is equal to twice the energy of the incident electron beam. 

Jets are reconstructed in the hadronic center-of-mass system 3 using the longitudinally boost 
invariant ^-algorithm and the E T -recombination scheme. The axis of each reconstructed 
jet is required to be within -1 < r\ = - ln(tan |) < 2.5 to ensure that the jets are well contained 
within the acceptance of the LAr calorimeter. Finally, a minimum transverse jet energy, E* T , of 
5 GeV is required. Demanding events with at least two jets which fulfill the criteria listed above 
yields a total sample of ~ 36 000 inclusive dijet events. 



4 Theoretical Predictions 

The available NLO QCD dijet and 3-jet programs provide the partonic final state of the hard 
subprocess, to which the chosen jet algorithm and selection can be applied. A variety of NLO 
dijet programs O0"ll3"Tl l3*2l have been shown to give comparable results 03113^1 l35l . Here we 
use a slightly modified version of DISENT OUIl in which the renormalization scale, // 2 , may be 

7 -. 2 

set to any linear combination of the two relevant scales, Q and E T , where the latter represents 
the mean transverse energy squared of the two hardest jets. The renormalization scale fij is set 

- 2 7 

to E T which, for most of the kinematic range under study, is larger than Q . The factorization 

- 2 

scale fi r is taken to be 70 GeV 2 , i.e. the average transverse jet energy squared, {E" T ), of the 
event sample 4 . The CTEQ6M (CTEQ6L) PDF parameterizations [36 1 are used for all NLO 
(LO) predictions shown. For NLO 3-jet production, the program NLOJET [031 is used. 

Theoretical predictions beyond the DGLAP collinear approach, which incorporate Iow-jc ef- 
fects by assuming different dynamics for the exchanged parton cascade, are available in Monte 
Carlo event generators. CCFM evolution, based on k t factorized unintegrated parton distri- 
butions, is implemented in the CASCADE generator [37 1 for initial state gluon showers. An 
alternative approach is provided by the ARIADNE Monte Carlo [38 1 program, which gener- 
ates non-fc r ordered parton cascades based on the color dipole model [39|. A LO Monte Carlo 
prediction, including effects due to the resolved hadronic structure of the virtual photon, and 
generating fc r ordered parton cascades as in the standard DGLAP approximation, is provided by 
RAPGAP RAPGAP can be run with ('direct+resolved') and without ('direct') a resolved 
photon contribution and the data are compared with both scenarios. RAPGAP ('direct') thus 
also allows a comparison with the standard DGLAP approach including full simulation of the 
hadronic final state. 

The LEPTO PT1 Monte Carlo program, which models only direct photon processes within 
the standard DGLAP approximation, and ARIADNE are used to estimate the hadronization 
corrections to be applied to the NLO predictions. All Monte Carlo models used here frag- 
ment the partonic final state according to the LUND string model [42 1 as implemented in JET- 
SET/PYTHIA El. 

3 Variables measured in the hadronic cms are marked by a '*'. 
4 DISENT does not allow to be varied with E* T event by event. 
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CASCADE ARIADNE RAPGAP LEPTO 



Version 


1 .U 


A 1 


Z.o 


O.J 


Proton PDF 


JS2001 El 
J2003 |46| 


CTEQ5L[47| 


CTEQ5L 


CTEQ5L 


Photon PDF 






SAS1D ED 




Renorm. scale yu 2 . 


* 2 2 

Pt +m q 


« 2 

Pt 


e 2 +V r 2 


Q 2 


Factor, scale ^u 2 


given by <mg. ordering 


. 2 

n 


n 2 +40* 2 




Underlying 
model 


CCFM 


Color 
dipole model 


DGLAP 
+ y-structure 


DGLAP 


Purpose 


Model comp. 


Model comp. 
QED/had. corr. 
detector corr. 


Model comp. 

QED corr. 
detector corr. 


had. corr. 



Table 1: Monte Carlo programs employed in the analysis. 



Higher order QED corrections are simulated using HERACLES |44|, which is directly in- 
terfaced to RAPGAP and via the DJANGO El program to ARIADNE. Both RAPGAP (direct) 
and ARIADNE are used to estimate the corrections for QED radiation and for detector effects 
as is outlined in the next section. 

The Monte Carlo programs used in the analysis are summarized in TableUJ including their 
basic settings. The background contribution from photoproduction events is estimated with the 
PHOJET [49 1 Monte Carlo program. 



5 Correction Procedure 

In order to compare data with theoretical predictions, the measured cross sections are corrected 
for detector acceptance and resolution, QED radiative effects and background contamination. 
In addition, hadronization corrections are applied to the NLO QCD calculations. The various 
correction factors are determined using the Monte Carlo models described above. These models 
reproduce the gross features of the jet data, as well as many characteristics of the final state, as 
shown in lESL However, none of the models gives a satisfactory description of all aspects 
of the hadronic final state. The most important discrepancies are found in the jet transverse 
momentum spectra. Differences between the Monte Carlo models are used to estimate the 
systematic uncertainties of this procedure. 

The corrections are applied to the data after statistical subtraction of the remaining pho- 
toproduction background. This contamination is estimated using the PHOJET Monte Carlo 
program. It is concentrated in the low-x region and is everywhere less than 4%. 

Detector and QED corrections are estimated using events generated with ARIADNE and 
RAPGAP (direct) and subjected to a full detector simulation and event reconstruction. The 
final correction factors are taken to be the average of the estimates from these two models. Half 
of the difference is included in the systematic uncertainty of the measurement. For the chosen 
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Source of error contributions Variation Uncertainty 



Experimental 


Hadronic energy scale 


4% 


7% 




SpaCal electromagnetic energy scale 


1% 


5% 




SpaCal hadronic energy scale 


7% 


2% 




Electron Polar angle measurement 


0.5 mrad 


2% 




Model uncertainty 




5 - 10% 




Photoproduction background 


30% 


1% 




Normalization uncertainty 




1.5% 


Theoretical 


Hadronization corrections 




5% 




Renormalization scale uncertainty 




10-30% 



Table 2: Summary of error contributions and the resulting typical uncertainties on the dijet 
cross section measurements (experimental) and the NLO predictions (theoretical). 



bins purities and stabilities are better than 40% for all data points. Here, the purity (stability) is 
defined as the number of dijet events which are both generated and reconstructed in a specific 
analysis bin, divided by the total number of dijet events that are reconstructed (generated) in 
that bin. The correction factors are in general between 0.8 and 1.2, but reach 1.8 at the lowest x 
and Q 2 values due to acceptance constraints in the backward calorimeter [ 28 1 . Additional minor 
corrections are applied to account for trigger inefficiencies. 

As mentioned before, hadronization corrections to the DISENT and NLOJET predictions 
are estimated using LEPTO and ARIADNE. The correction factors are determined by compar- 
ing the cross sections calculated from the hadronic final state (hadron level) with those predicted 
from the partonic final state (LO and QCD parton showers) prior to the hadronization step. They 
are obtained by taking the average of the estimates derived from LEPTO and from ARIADNE. 
When applied to the NLO predictions, these corrections allow for comparisons between data 
and theory at the hadron level. The correction factors lower the NLO predictions by typically 
10%. Half of the difference between the two models is taken as the systematic error on the 
hadronization correction. 



6 Systematic Uncertainties 

The different error sources and the corresponding uncertainties on the dijet cross section mea- 
surements are summarized in Table EJ The theoretical uncertainties on the NLO predictions, 
given by the errors on the hadronization corrections and the renormalization scale uncertainty, 

9 - 2 -.2 

are also listed. The latter is estimated by varying yr r between E T /4 and 4 E T . 

One of the most important error contributions arises from the uncertainty in the hadronic 
energy measurement used in the jet reconstruction. This scale uncertainty was estimated to be 
4% and leads to an uncertainty of typically 7% on the dijet cross section measurement, with 
values increasing up to 20% at large transverse jet energies. The uncertainty of the electromag- 
netic energy scale of the SpaCal is 1% and leads to an error on the dijet cross sections of 5% in 
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most of the phase space, reaching ~10% at large x, where in some bins it constitutes the largest 
contribution to the total systematic error. The influence of the hadronic energy scale uncertainty 
of the SpaCal of 7% is of minor importance. It only enters in the determination of ^(it, - p z j) 
and gives a 2% contribution to the final measurement error. An error of similar size arises from 
the polar angle measurement of the scattered electron. 

The differences between the correction factors when using different Monte Carlo models 
lead to an error contribution of ~ 5 to 10% throughout the analyzed phase space. The 30% 
uncertainty on the absolute normalization of the yp-background contributes up to 1% to the 
systematic error on the dijet cross sections. 

The total systematic error is determined by summing the individual contributions in quadra- 
ture. A 1.5% normalization uncertainty due to the luminosity measurement is not included in 
the quoted systematic errors on the cross sections presented here. 



7 Results 

7.1 Inclusive Dijet Cross Sections 

All measured dijet cross sections are presented after correcting for detector and radiative effects. 
They are given multi-differentially as a function of x, Q 2 and several dijet observables and 
are compared with the DISENT NLO calculations after applying hadronization corrections. 
NLO calculations of dijet observables become sensitive to soft gluon radiation when symmetric 
selection criteria on the transverse jet energies are applied 1 3"31 15011511 . Thus, in addition to 
the requirement E* T > 5 GeV for the two highest transverse momentum jets, an additional 
requirement on the most energetic jet, E* T x = E* Tmax > (5 + A) GeV, is necessary. This avoids 
regions of phase space in which NLO predictions become unreliable. Figure El shows the dijet 
cross section as a function of the parameter A in bins of x and Q 2 . Within the theoretical 
uncertainties, good agreement between the data and the NLO predictions is found for all values 
of A. However, an unphysical reduction in the NLO calculation occurs for A < 1 GeV. For 
comparison the figure also presents the LO DISENT prediction at the parton level. The large 
differences between the LO and the NLO predictions, as well as the large scale uncertainties in 
the NLO predictions, indicate the need for higher order contributions, especially at low x and 
low Q 2 . 

Figure |3] shows the dijet cross section as a function of Bjorken-x in intervals of Q 2 for 
fixed A = 2 GeV. The data show a significant increase towards low x, which is consistent 
with the strong rise of the gluon density observed in low-x structure function measurements at 
HERA 1261521 . No deviation of the data from calculations using the conventional NLO DGLAP 
approach is found. 5 The scale uncertainties are sizable and increase towards low x. 

Measurements of the dijet cross section for A = 2 GeV as a function of E* Tmax and \Arj*\ 
in bins of Bjorken-x and Q 2 are shown in Figures HI and Within the quoted uncertainties, 

5 Minor differences at low x and Q 2 as reported in 1 19 1 are still observed when using CTEQ4M |53|, an older 
parameterization of the parton distribution functions. 
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good agreement between data and NLO calculations is observed, even at small values of £!„ 
and |A?7*|, i.e. in a kinematic region in which effects due to low-x dynamics should be most 
prominent ITP71 . The level of agreement is more visible in Figures^ and |7J which show the ratio 
of the data to the NLO predictions as a function of E* T max and | A77* | in bins of Bjorken-x and Q 2 . 

The measured dijet cross sections are summarized in Tables|3]to|5] in addition, all dijet cross 
sections shown have been normalized to the total inclusive cross section for each bin. These 
dijet rates, R 2 = A^ di j et /A^ DI s, are listed in Tables|6]to|E] 

7.2 Azimuthal Jet Separation 

Insight into low-x dynamics can be gained from inclusive dijet data by studying the behavior of 
events with a small azimuthal separation, A(p* , between the two hardest jets as measured in the 
hadronic center-of-mass system H181ll7lfT6l . Partons entering the hard scattering process with 
negligible transverse momentum, k t , as assumed in the DGLAP formalism, produce at leading 
order a back-to-back configuration of the two outgoing jets with A(p* ~ 180°. Azimuthal jet 
separations different from 180° occur due to higher order QCD effects. However, in models 
which predict a significant proportion of partons entering the hard process with large k t , the 
number of events with small A(f>* increases. This is the case for the BFKL and CCFM evolution 
schemes. As an illustration, Figure[S]shows the uncorrected A(p* distribution for several intervals 
in Q 2 . The expected steeply falling spectrum is observed with a tail extending to small values 
of A0*. This behavior is broadly reproduced by the Monte Carlo programs RAPGAP (direct) 
and ARIADNE, which are used to correct the data and estimate model uncertainties in the same 
manner as for the differential cross section measurements. 

Large migrations connected with the limited hadronic energy resolution make an extraction 
of the dijet cross section at small A</>* rather difficult. Thus the ratio 

f° N di]et (Act>*,x, Q 2 )dAcf>* 

S ( a ) = TTro^ ' 

X jV dije t(A0*,x,<2 2 )dA0* 

of the number of events A^et with an azimuthal jet separation of A(f>* < a relative to all dijet 
events is measured, as proposed in ifTol . This variable is also directly sensitive to Iow-jc effects. 
For the analysis presented here a = 120° is chosen, which results in a purity of around 45% 
independently of x and Q 2 and systematic uncertainties of similar size to those on the cross 
section measurements. 

Figure presents the S distribution for a = 120° as a function of x for different Q 2 . The 
measured S values are summarized in Tabled For the chosen a, the measured values of S are 
of the order of 5% and increase with decreasing x. This rise of S is most prominent in the lowest 
Q 2 bin, where the lowest values of x are reached. The NLO dijet QCD calculations predict S 
values of only ~1% and show no rise towards low x. Low values of S are expected for the 
NLO dijet predictions, since without any restrictions in acceptance, the two most energetic jets 
should always be separated by more than Acp* = 120°. However, since selection criteria have 
to be applied to match the experimental conditions, non-vanishing S values arise, due to event 
topologies for which some of the jets lie outside the analyzed phase space. In the same figure 
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NLO 3-jet predictions are also shown. These give a good description of the data at large Q 2 and 
large x, but still fail to describe the increase towards low x, particularly in the lowest Q 2 range. 

An informative comparison of the measured S distribution with theory is also provided 
by models with different implementations of higher order QCD effects. Within the DGLAP 
approach such a model is provided for example by RAPGAP (direct). As shown in Figure fTOl 
RAPGAP (direct) predicts a much larger ratio S than the NLO dijet calculation. However, it 
still fails to describe the data in the Iow-jc, low-<2 2 region. An improved description is achieved 
when resolved photon processes are included in RAPGAP (direct+resolved). Even with direct 
and resolved photon contributions included, RAPGAP fails to describe the data at very low x 
and Q 2 . Note that to obtain this overall level of agreement, it was necessary to choose a rather 

large scale, i.e. p 2 = Q 2 + AE' T , in order to get a large enough resolved photon contribution. 

If the observed discrepancies are due to the influence of non-fc r ordered parton emissions, 
models based on the color dipole model or CCFM evolution may provide a better description of 
the ratio S . In Figure [TT]the data are therefore compared with the predictions of the ARIADNE 
and CASCADE Monte Carlo programs. For CASCADE, the two predictions presented are 
based on the JS2001 [37 1 and set 2 of the J2003 [46 1 unintegrated parton distributions, which 
differ in the way the small k, region is treated. For J2003 set 2, the full splitting function, 
including the non-singular term, is used, in contrast to JS2001 , for which only the singular terms 
were considered. Whereas the prediction for S using JS2001 lies significantly above the data, 
that based on set 2 of J2003 describes the data rather well. Note that both PDFs describe the HI 
structure function data Hence, the measurement of the ratio S improves the sensitivity to 
the details of the unintegrated gluon distribution. A good description of the S distribution at low 
x and low Q 2 is also provided by the color dipole model incorporated in ARIADNE. However, 
at higher Q 2 the ARIADNE prediction falls below the measured S values. 

8 Conclusion 

Inclusive dijet production in deep inelastic ep scattering is measured in the kinematic range 
5 < Q 2 < 100 GeV 2 , 10" 4 < x < 10" 2 and 0.1 < y < 0.7. Multi-differential cross section data 
are compared with NLO QCD predictions and no significant deviations are observed within the 
experimental and theoretical uncertainties. In the kinematic range studied, the next-to-leading 
order DGLAP approach thus provides an adequate theory for predicting ep dijet cross sections 
as a function of Bjorken-x, Q 2 , Ej max and |A?f|. 

NLO dijet QCD calculations predict values that are much too low for the ratio, S , of events 
with a small azimuthal separation of the two highest transverse momentum jets with respect 
to the total number of inclusive dijet events. The additional hard emission, provided by the 
NLO 3-jet calculations, considerably improves the description of the data, but is insufficient at 
low x and low Q 2 . A similar description of the data is provided by RAPGAP, a DGLAP-based 
QCD model, which matches LO matrix elements for direct and resolved processes to fc r ordered 
parton cascades. A good description of the measured ratio S at low x and Q 2 is given by the 
ARIADNE program, which generates non-^-ordered parton cascades using the color dipole 
model. Predictions based on the CCFM evolution equations and k, factorized unintegrated 
gluon densities are provided by the CASCADE Monte Carlo program. Large differences are 
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found between the predictions for two different choices of the unintegrated gluon density, both 
of which describe the HI structure function data and one of which gives a good description of 
S . This measurement thus provides a significant constraint on the unintegrated gluon density. 
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Figure 1 : Leading order diagrams for dijet production in ep scattering, (a) photon-gluon fusion 
and (b) QCD-Compton process, (c) parton cascade diagram: k, denotes the transverse momenta 
of the exchanged gluons, x g the fractional longitudinal momentum of the gluon taking part in 
the hard process and x is the Bjorken scaling variable. 
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Figure 2: Inclusive dijet cross section as given in Tabled multiplied by (x) and averaged over x 
and Q 2 , as a function of A, defined by the requirement E* Tl = E* T max > (5 + A) GeV. Here (x) and 
(Q 2 ) are the mean values of Bjorken-jc and Q 2 for fully inclusive events in a given bin. The data 
are shown together with their statistical uncertainties (inner error bars) and their statistical and 
systematic uncertainties added in quadrature (outer error bars). They are compared with NLO 
(LO) dijet QCD predictions using the CTEQ6M (CTEQ6L) parton distribution functions. The 
NLO predictions are corrected for hadronization effects. The outer light error band includes the 
quadratic sum of hadronization (dark error band) and renormalization scale uncertainties on the 
NLO predictions. 
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Figure 3: Inclusive dijet cross section averaged over Q 2 and Bjorken-x for A = 2 GeV (see text). 
The data are plotted at the center of each bin and are shown with their statistical uncertainties 
(inner error bars) and their statistical and systematic uncertainties added in quadrature (outer 
error bars). They are compared with NLO (LO) dijet QCD predictions using the CTEQ6M 
(CTEQ6L) parton distribution functions. The NLO predictions are corrected for hadronization 
effects. The outer light error band includes the quadratic sum of hadronization (dark error band) 
and renormalization scale uncertainties on the NLO predictions. 
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Figure 4: Inclusive dijet cross section for A = 2 GeV averaged over Bjorken-jc, Q 2 and 
as given in Table El compared with NLO dijet QCD predictions using the CTEQ6M parton 
distribution functions. The data are plotted at the center of each bin and are shown together with 
their statistical uncertainties (inner error bars) and their statistical and systematic uncertainties 
added in quadrature (outer error bars). The NLO predictions are corrected for hadronization 
effects. The outer light error band includes the quadratic sum of hadronization (dark error 
band) and renormalization scale uncertainties on the NLO predictions. 
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Figure 5: Inclusive dijet cross section for A = 2 GeV averaged over Bjorken-x, Q 2 and the 
pseudorapidity distance \Atj*\ between the dijets as given in Table compared with NLO dijet 
QCD predictions using the CTEQ6M parton distribution functions. The data are plotted at the 
center of each bin and are shown together with their statistical uncertainties (inner error bars) 
and their statistical and systematic uncertainties added in quadrature (outer error bars). The 
NLO predictions are corrected for hadronization effects. The outer light error band includes the 
quadratic sum of hadronization (dark error band) and renormalization scale uncertainties on the 
NLO predictions. 
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Figure 6: The ratio of the measured inclusive dijet cross section for A = 2 GeV to the theo- 
retical prediction in bins of Bjorken-x, Q 2 and E* Tmax . The data are shown together with their 
statistical uncertainties (inner error bars) and their statistical and systematic uncertainties added 
in quadrature (outer error bars). They are compared with NLO dijet QCD predictions using the 
CTEQ6M parton distribution functions. The theoretical errors are given by the light error band 
representing the quadratic sum of the hadronization (dark error band) and renormalization scale 
uncertainties. 
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Figure 7: The ratio of the measured inclusive dijet cross section for A = 2 GeV to the the- 
oretical prediction in bins of Bjorken-x, Q 2 and |Aj7*|. The data are shown together with their 
statistical uncertainties (inner error bars) and their statistical and systematic uncertainties added 
in quadrature (outer error bars). They are compared with NLO dijet QCD predictions using the 
CTEQ6M parton distribution functions. The theoretical errors are given by the light error band 
representing the quadratic sum of the hadronization (dark error band) and renormalization scale 
uncertainties. 
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Figure 8: Normalized measured event distributions as a function of the azimuthal separation, 
A(f>*, between the two highest transverse momentum jets. The data are shown for different 
regions in Q 2 and are compared with two Monte Carlo models, ARIADNE and RAPGAP. 
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Figure 9: Ratio S of the number of events with a small azimuthal jet separation (Acf>* < 120°) 
between the two highest transverse momentum jets with respect to the total number of inclusive 
dijet events, as a function of Bjorken-jc and Q 2 . The data are plotted at the center of each bin 
and are shown together with their statistical uncertainties (inner error bars) and their statistical 
and systematic uncertainties added in quadrature (outer error bars). The data are compared with 
NLO QCD predictions for dijet and 3-jet production using the CTEQ6M parton distribution 
functions. The theoretical errors are given by the light error band representing the quadratic 
sum of the hadronization (dark error band) and renormalization scale uncertainties. 
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Figure 10: Ratio S of the number of events with a small azimuthal jet separation (A0* < 120°) 
between the two highest transverse momentum jets with respect to the total number of inclusive 
dijet events, as a function of Bjorken-jc and Q 2 . The data are plotted at the center of each bin 
and are shown together with their statistical uncertainties (inner error bars) and their statistical 
and systematic uncertainties added in quadrature (outer error bars). The data are compared with 
predictions from the RAPGAP generator, both with direct photons alone (full line) and with 
direct and resolved contributions (dashed line). 
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Figure 11: Ratio S of the number of events with a small azimuthal jet separation (A0* < 120°) 
between the two highest transverse momentum jets with respect to the total number of inclusive 
dijet events, as function of Bjorken-* and Q 2 . The data are plotted at the center of each bin 
and are shown together with their statistical uncertainties (inner error bars) and their statistical 
and systematic uncertainties added in quadrature (outer error bars). They are compared with 
predictions from a model based on CCFM evolution (CASCADE) and using two different k t - 
unintegrated gluon distribution functions, JS2001 (full line) and set 2 of J2003 (dashed line). In 
addition, the data are compared with predictions based on the color dipole model (ARIADNE), 
which produces £ r unordered parton showers. 
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Q 2 

[GeV 2 ] 


x/10~ 4 


A 
[GeV] 


[pb/GeV 2 ] 


$stat &syst. 

[%] [%] 


5.0 - 10.0 


1.0 - 1.7 



1 
2 
4 
7 


6.1 • 10 5 
5.6 ■ 10 5 

4.6 • 10 5 

2.7 ■ 10 5 

1.2 • 10 5 


3 8 
3 7 

3 7 

4 8 
6 9 




1.7 - 3.0 



1 
2 
4 
7 


4.8 ■ 10 3 

4.3 • 10 5 

3.4 ■ 10 5 
1.8 • 10 5 
6.6 ■ 10 4 


3 9 
3 10 

3 9 

4 9 

5 11 




3.0 - 5.0 



1 
2 
4 
7 


2.3 • 10 3 
2.0 ■ 10 5 
1.5 ■ 10 5 
8.0 ■ 10 4 
2.9 • 10 4 


4 12 
4 12 

4 11 

5 11 
8 13 




5.0 - 10.0 



1 
2 
4 
7 


5.2 • 10 4 
4.6 ■ 10 4 
3.4 ■ 10 4 
1.6 ■ 10 4 
5.4 • 10 3 


5 14 
5 13 
5 13 
7 12 
10 15 


10.0 - 15.0 


1.7 - 3.0 



1 
2 
4 
7 


1.7- 10 s 
1.6- 10 5 
1.4- 10 5 
8.7 ■ 10 4 
3.7 • 10 4 


4 8 
4 8 

4 9 

5 9 
7 8 




3.0 - 5.0 



1 
2 
4 
7 


1.3 • 10 3 
1.2- 10 5 
9.1 ■ 10 4 
4.8 • 10 4 
2.0 • 10 4 


3 10 

3 9 

4 9 

5 10 
7 11 




5.0 - 10.0 



1 
2 
4 
7 


5.9 • 10 4 

5.2 ■ 10 4 
4.0 • 10 4 
2.0 • 10 4 

7.3 • 10 3 


3 13 

3 12 

4 12 

5 12 
7 13 




10.0 - 18.0 



1 
2 
4 
7 


1.1 • 10 4 
8.8 ■ 10 3 
6.4 • 10 3 
3.3 ■ 10 3 
9.3 • 10 2 


6 12 

7 13 
7 15 
9 15 
17 22 



Table 3: Inclusive dijet cross section averaged over the regions indicated in x and Q 2 , for dif- 
ferent values of A as shown in Figures|2]and|3] The measurement is restricted to values of the 
inelasticity variable y between 0.1 < y < 0.7 and to values of the polar angle of the scattered 
electron between 156° < 9 < 175°. 
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Q 2 
[GeV 2 ] 


x/10- 4 


A 
[GeV] 


[pb/GeV 2 ] 


Sstat Ssyst. 
[%] [%] 


15.0 - 20.0 


3.0 - 5.0 



1 
2 
4 
7 


7.8- 10 4 

7.6 ■ 10 4 

6.4 ■ 10 4 

3.7 ■ 10 4 

1.5 - 10 4 


6 9 
4 9 

4 9 

5 9 

7 10 




5.0 - 10.0 



1 
2 
4 
7 


4.5 • 10 4 

4.0 • 10 4 

3.1 • 10 4 

1.6 ■ 10 4 
5.8 • 10 3 


4 11 
4 12 

4 11 

5 11 
8 12 




10.0 - 22.0 



1 
2 
4 
7 


1.1 • 10 4 
9.9 ■ 10 3 
7.6 ■ 10 3 
3.8 • 10 3 
1.3 • 10 3 


5 12 
5 13 
5 12 
7 13 
11 15 


20.0 - 30.0 


3.0 - 5.0 



1 
2 
4 
7 


2.5 • 10 4 
2.4 • 10 4 
2.1 ■ 10 4 
1.4 ■ 10 4 
6.7 • 10 3 


5 8 
5 8 

5 9 

6 8 
9 9 




5.0 - 10.0 



1 
2 
4 
7 


2.9 • 10 4 
2.7 • 10 4 
2.1 ■ 10 4 
1.2- 10 4 
5.1 • 10 3 


3 9 
3 9 

3 10 

4 10 
6 10 




10.0 - 33.0 



1 
2 
4 
7 


6.2 • 10 3 
5.5 ■ 10 3 

4.1 ■ 10 3 

2.2 • 10 3 

9.3 • 10 2 


3 11 
3 11 

3 11 

4 12 
7 13 


30.0 - 50.0 


5.0 - 10.0 



1 
2 
4 
7 


1.0 • 10 4 

9.5 ■ 10 3 
8.2- 10 3 
5.4- 10 3 

2.6 • 10 3 


3 8 

3 8 

4 10 
4 9 
6 10 




10.0 - 55.0 



1 
2 
4 
7 


2.8 • 10 3 
2.5 • 10 3 

2.0 ■ 10 3 

1.1 ■ 10 3 
4.7 • 10 2 


2 10 

2 10 

3 11 
3 12 
5 13 


50.0 - 100.0 


10.0 - 25.0 



1 
2 
4 
7 


1.4- 10 3 
1.3 ■ 10 3 
1.1 ■ 10 3 
6.5 ■ 10 2 
2.8- 10 2 


3 9 
3 9 

3 9 

4 10 
6 13 




25.0 - 100.0 



1 
2 
4 
7 


4.5 • 10 2 

4.0 • 10 2 
3.2- 10 2 
1.8 ■ 10 2 

7.1 • 10 1 


3 9 
3 10 

3 11 

4 11 
6 17 



Table |3] continued . 
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Q 2 

[GeV 2 ] 


x/10- 4 


T,max 
[GeV] 


d 3 cr 
dQ 2 dxdE* 

[pb/GeV 3 ] 


Sstat Ssyst. 
[%] [%] 


5.0 - 10.0 


1.0 - 1.7 


7.0 - 12.0 
12.0 - 20.0 
20.0 - 30.0 
30.0 - 60.0 


6.8- 10 4 
1.2- 10 4 
1.3 ■ 10 3 

H ^> 1 Al 

7.2 • 10 


4 9 
6 12 
13 16 

11 11 




1.7 - 3.0 


7.0 - 12.0 
12.0 - 20.0 

1A A OA A 

20.0 - 30.0 
30.0 - 60.0 


5.5 • 10 4 
6.8- 10 3 

£. C\ 1 a2 

6.9 • 10 
3.0- 10 1 


3 11 
5 14 
13 17 
35 25 




3.0 - 5.0 


7.0 - 12.0 

1 i a a 

12.0 - 20.0 
20.0 - 30.0 


2.5 ■ 10 4 

i 1 1 a3 

3.1 ■ 10 
2.9- 10 2 


4 12 

O 1 "7 

8 17 
21 21 




5.0 - 10.0 


7.0 - 12.0 
12.0 - 20.0 

OA A 1 A A 

20.0 - 30.0 


5.6- 10 3 
6.0- 10 2 

1 A 1 Al 

3.0 • 10 


5 15 

10 17 

11 1A 

33 30 


10.0 - 30.0 


1.7 - 3.0 


7.0 - 12.0 
12.0 - 20.0 
20.0 - 30.0 

OA A /" /\ r\ 

30.0 - 60.0 


5.5- 10 3 
1.1 • 10 3 
1.0- 10 2 

C A 1 AO 

5.4 • 10 


6 11 

7 12 
18 14 

36 22 




3.0 - 5.0 


7.0 - 12.0 
12.0 - 20.0 
20.0 - 30.0 

OA A /" /\ r\ 

30.0 - 60.0 


7.4- 10 3 
1.2- 10 3 

1.5- 10 2 

1 A 1 Al 

1.0 • 10 


3 11 
5 13 
11 18 
20 25 




5.0 - 1.0 


7.0 - 12.0 
12.0 - 20.0 
20.0 - 30.0 

1 A A H A A 

30.0 - 60.0 


4.6 ■ 10 3 
6.2- 10 2 
6.1 ■ 10 1 
1.4 ■ 10° 


2 12 
4 15 
11 18 
36 24 




1.0 - 33.0 


7.0 - 12.0 
12.0 - 20.0 
20.0 - 30.0 

■1 A A £ A A 

30.0 - 60.0 


5.7- 10 2 

7.8- 10 1 
4.6 ■ 10° 

A 1 1 AO 

0.1 • 10 


3 14 
6 16 
18 26 
69 23 


30.0 - 100.0 


5.0 - 10.0 


7.0 - 12.0 
12.0 - 20.0 
20.0 - 30.0 

OA A £A A 

30.0 - 60.0 


3.3- 10 2 
7.8- 10 1 
9.7 • 10° 

AO 1 AO 

0.3 • 10 


4 11 
6 14 
17 24 
39 32 




10.0 - 100.0 


7.0 - 12.0 
12.0 - 20.0 
20.0 - 30.0 
30.0 - 60.0 


9.3- 10 1 
1.5- 10 1 
1.3- 10° 
5.0- 10~ 2 


2 11 

3 17 
9 18 
26 30 



Table 4: Inclusive dijet cross section averaged over the regions indicated in x, Q 2 and E* Tmax 
for A = 2 GeV as shown in FigureHJ The measurement is restricted to values of the inelastic- 
ity variable y between 0.1 < y < 0.7 and to values of the polar angle of the scattered electron 
between 156° < 6 < 175°. 
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Q 2 

[Uev J 


x/10- 4 


|A7fl 


d Q 2 dxd\ A?]*\ 
[pb/uev J 


Ostat $syst. 

l%J l%J 


5.0 - 10.0 


1.0 - 1.7 


0.0 - 0.7 
0.7 - 1.4 
1.4 - 2.1 
2.1 - 2.8 
2.8 - 3.5 


2.3 • 10 5 
1.9 10 5 
1.3 ■ 10 5 
7.0 • 10 4 
2.0 • 10 4 


4 10 

5 11 

6 9 
9 9 
18 13 




1.7 - 3.0 


0.0 - 0.7 
0.7 - 1.4 
1.4 - 2.1 
2.1 - 2.8 
2.8 - 3.5 


1.7 • 10 5 
1.4- 10 5 

1.0 ■ 10 s 

5.1 ■ 10 4 
1.1 ■ 10 4 


4 12 
4 12 
6 14 
8 12 
17 30 




3.0 - 5.0 


0.0 - 0.7 
0.7 - 1.4 
1.4 - 2.1 
2.1 - 2.8 


8.5 • 10 4 
6.9 ■ 10 4 

4.6 ■ 10 4 

1.7 ■ 10 4 


5 14 

6 15 
8 13 
13 17 




5.0 - 10.0 


0.0 - 0.7 
0.7 - 1.4 
1.4 - 2.1 


2.0 • 10 4 
1.7 ■ 10 4 
8.6 ■ 10 3 


7 14 

8 16 
12 22 


10.0 - 30.0 


1.7 - 3.0 


0.0 - 0.7 
0.7 - 1.4 
1.4 - 2.1 
2.1 - 2.8 
2.8 - 3.5 


2.0 • 10 4 
1.5- 10 4 
1.2- 10 4 
5.0- 10 3 
1.8 • 10 3 


5 13 

6 10 

7 10 
11 13 
20 31 




3.0 - 5.0 


0.0 - 0.7 
0.7 - 1.4 
1.4 - 2.1 
2.1 - 2.8 
2.8 - 3.5 


2.5 • 10 4 
1.9 ■ 10 4 
1.4- 10 4 
9.4 ■ 10 3 
2.0 • 10 3 


4 12 

4 11 

5 11 
7 11 
13 13 




5.0 - 1.0 


0.0 - 0.7 
0.7 - 1.4 
1.4 - 2.1 
2.1 - 2.8 
2.8 - 3.5 


1.5- 10 4 
1.2- 10 4 
8.1 ■ 10 3 
4.0 ■ 10 3 
8.8 • 10 2 


3 12 
3 12 
5 15 
7 14 
16 30 




1.0 - 33.0 


0.0 - 0.7 
0.7 - 1.4 
1.4 - 2.1 
2.1 - 2.8 


2.1 • 10 3 
1.8 ■ 10 3 
9.5 ■ 10 2 

2.2 • 10 2 


4 13 
4 15 
7 16 
13 26 


30.0 - 100.0 


5.0 - 10.0 


0.0 - 0.7 
0.7 - 1.4 
1.4 - 2.1 
2.1 - 2.8 
2.8 - 3.5 


1.1 • 10 3 
1.1 ■ 10 3 
6.8 • 10 2 
3.4 • 10 2 
1.6 ■ 10 2 


6 11 

6 14 

7 12 
11 13 
21 24 




10.0 - 100.0 


0.0 - 0.7 
0.7 - 1.4 
1.4 - 2.1 
2.1 - 2.8 
2.8 - 3.5 


3.4 • 10 2 

2.7 ■ 10 2 

1.8 ■ 10 2 
6.1 ■ 10 1 
9.1 • 10° 


3 14 

3 11 

4 13 
6 14 
15 21 



Table 5: Inclusive dijet cross section averaged over the regions indicated in x, Q 2 and |A?7*| 
for A = 2 GeV as shown in Figure|5] The measurement is restricted to values of the inelastic- 
ity variable y between 0.1 < y < 0.7 and to values of the polar angle of the scattered electron 
between 156° < 9 < 175°. 
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Q 2 
[GeV 2 ] 


x/10- 4 


A 
[GeV] 


^ 2 - N DIS 


$stat &syst. 

[%] [%] 


5.0 - 10.0 


1.0 - 1.7 



1 
2 
4 
7 


0.031 
0.028 
0.023 
0.014 
0.006 


4 9 
4 9 

3 9 

4 9 
6 10 




1.7 - 3.0 



1 
2 
4 
7 


0.027 
0.024 
0.019 
0.010 
0.004 


3 9 
3 9 

3 9 

4 10 

5 12 




3.0 - 5.0 



1 
2 
4 
7 


0.021 
0.019 
0.014 
0.007 
0.003 


4 12 
4 12 

4 11 

5 11 
8 13 




5.0 - 10.0 



1 
2 
4 
7 


0.015 
0.013 
0.010 
0.004 
0.002 


5 15 
5 14 
5 14 
7 14 
10 18 


10.0 - 15.0 


1.7 - 3.0 



1 
2 
4 
7 


0.040 
0.037 
0.031 
0.020 
0.009 


4 8 
4 8 

4 9 

5 10 
7 9 




3.0 - 5.0 



1 
2 
4 
7 


0.035 
0.031 
0.024 
0.013 
0.005 


3 10 

3 9 

4 9 

5 10 
7 11 




5.0 - 10.0 



1 
2 
4 
7 


0.027 
0.024 
0.018 
0.009 
0.003 


3 13 

4 12 

4 13 

5 13 
7 14 




10.0 - 18.0 



1 
2 
4 
7 


0.019 
0.016 
0.011 
0.006 
0.002 


6 15 

7 14 
7 18 
9 20 
17 24 



Table 6: Dijet rate R 2 averaged over the regions indicated in Q 2 and x for different values of A. 
The measurement is restricted to values of the inelasticity variable y between 0.1 < y < 0.7 and 
to values of the polar angle of the scattered electron between 156° < 9 < 175°. 
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Q 2 
[GeV 2 ] 


x/10- 4 


A 
[GeV] 


r> jet 


Sstat Ssyst. 
[%] [%] 


15.0 - 20.0 


3.0 - 5.0 



1 
2 
4 
7 


0.043 
0.042 
0.035 
0.020 
0.008 


6 9 
4 8 

4 9 

5 9 

7 10 




5.0 - 10.0 



1 
2 
4 
7 


0.041 
0.036 
0.028 
0.015 
0.005 


4 10 
4 12 

4 11 

5 10 
8 12 




10.0 - 22.0 



1 
2 
4 
7 


0.027 
0.024 
0.019 
0.009 
0.003 


5 13 
5 13 
5 12 
7 14 
11 17 


20.0 - 30.0 


3.0 - 5.0 



1 
2 
4 
7 


0.056 
0.052 
0.046 
0.030 
0.015 


5 8 

5 8 

6 8 

7 9 
9 10 




5.0 - 10.0 



1 
2 
4 
7 


0.054 
0.049 
0.040 
0.023 
0.010 


3 9 
3 9 

3 9 

4 10 
6 10 




10.0 - 33.0 



1 
2 
4 
7 


0.038 
0.033 
0.025 
0.014 
0.006 


3 11 

3 11 

4 12 
4 13 
7 15 


30.0 - 50.0 


5.0 - 10.0 



1 
2 
4 
7 


0.068 
0.064 
0.056 
0.036 
0.018 


3 8 

3 8 

4 9 
4 9 
6 10 




10.0 - 55.0 



1 
2 
4 
7 


0.056 
0.051 
0.041 
0.023 
0.009 


2 10 

2 10 

3 11 
3 12 
5 13 


50.0 - 100.0 


10.0 - 25.0 



1 
2 
4 
7 


0.087 
0.080 
0.067 
0.041 
0.018 


3 9 
3 8 

3 9 

4 10 
6 13 




25.0 - 100.0 



1 
2 
4 
7 


0.072 
0.065 
0.052 
0.029 
0.011 


3 9 
3 10 

3 11 

4 12 
6 18 



Table |6] continued . 



Q 2 

[GeV 2 ] 


x/10~ 4 


[GeV] 


dR 2 
dE* 

l ,max 

[GeV- 1 ] 


Sstat &syst. 
[%] [%] 


5.0 - 10.0 


1.0 - 1.7 


7.0 - 12.0 
12.0 - 20.0 
20.0 - 30.0 
30.0 - 60.0 


3.5 • 10~ 3 
4.9 • 10~ 4 
6.8 • 10~ 5 
3.6- 10~ 6 


3 10 
6 14 
13 18 
27 27 




1.7 - 3.0 


7.0 - 12.0 
12.0 - 20.0 
20.0 - 30.0 
30.0 - 60.0 


3.1 • 10~ 3 
3.1 • 10~ 4 
4.0 • 10~ 5 
1.7- 10~ 6 


3 11 
5 14 
13 18 
35 25 




3.0 - 5.0 


7.0 - 12.0 
12.0 - 20.0 
20.0 - 30.0 


2.3 • 10~ 3 
2.3 • 10~ 4 
2.6 • 10~ 5 


4 13 
8 17 
21 22 




5.0 - 10.0 


7.0 - 12.0 
12.0 - 20.0 
20.0 - 30.0 


1.6- 10~ 3 
1.4- 10~ 4 
7.8 • 10~ 6 


5 16 
10 20 
34 36 


10.0 - 30.0 


1.7 - 3.0 


7.0 - 12.0 
12.0 - 20.0 
20.0 - 30.0 
30.0 - 60.0 


4.6 • 10~ 3 
7.2 • 10~ 4 
8.6- 10~ 5 
4.5 • 10~ 6 


5 11 
7 13 
18 15 
36 23 




3.0 - 5.0 


7.0 - 12.0 
12.0 - 20.0 
20.0 - 30.0 
30.0 - 60.0 


4.6 • 10~ 3 
5.9 • 10~ 4 
9.2 • 10~ 5 
6.4 • 10~ 6 


3 11 
5 13 
11 18 
20 25 




5.0 - 1.0 


7.0 - 12.0 
12.0 - 20.0 
20.0 - 30.0 
30.0 - 60.0 


4.2 • 10~ 3 
4.6 • 10~ 4 
5.6- 10- 5 
1.3- 10~ 6 


2 12 
4 15 
11 19 
36 25 




1.0 - 33.0 


7.0 - 12.0 
12.0 - 20.0 
20.0 - 30.0 
30.0 - 60.0 


3.1 • 10~ 3 

3.4 • 10~ 4 

2.5 • 10~ 5 
7.0- 10~ 7 


3 14 
6 17 
18 27 
69 25 


30.0 - 100.0 


5.0 - 10.0 


7.0 - 12.0 
12.0 - 20.0 
20.0 - 30.0 
30.0 - 60.0 


7.7 • 10~ 3 
1.5 - 10- 3 
2.3 • 10~ 4 
7.3 • 10~ 6 


4 11 
6 13 
17 23 
39 32 




10.0 - 100.0 


7.0 - 12.0 
12.0 - 20.0 
20.0 - 30.0 
30.0 - 60.0 


7.4 • 10~ 3 

9.5 • 10~ 4 
1.1 • 10~ 4 
4.0 • 10~ 6 


2 11 

3 17 
9 19 
26 31 



Table 7: Dijet rate R 2 averaged over the regions indicated in x, Q 2 and E* Tmax for A = 2 GeV. 
The measurement is restricted to values of the inelasticity variable y between 0.1 < y < 0.7 and 
to values of the polar angle of the scattered electron between 156° < 6 < 175°. 
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O 2 




x/10~ 4 


IAn*l 




dR 2 
d\Ar,*\ 


'J stat 


u syst. 




[GeV 2 ] 












[%] 


[%] 


5.0 


10.0 


1.0 


1.7 


0.0 - 


0.7 


1.2- 10~ 2 


5 


12 










0.7 - 


1.4 


9.5 • 10~ 3 


5 


12 










1.4 - 


2.1 


6.5 • 10~ 3 


7 


11 










2.1 - 


2.8 


3.5 • 10~ 3 


9 


10 










2.8 - 


3.5 


9.9 • 10~ 4 


18 


13 






1.7 


3.0 


0.0 - 


0.7 


9.6 • 10 3 


4 


12 










0.7 - 


1.4 


8.0- 10~ 3 


5 


12 










1.4 - 


2.1 


5.8 • 10~ 3 


6 


13 










2.1 - 


2.8 


2.9 • 10~ 3 


8 


12 










2.8 - 


3.5 


6.3 • 10- 4 


17 


31 






3.0 


5.0 


0.0 - 


0.7 


7.7 • 10~ 3 


5 


14 










0.7 - 


1.4 


6.3 • 10~ 3 


6 


15 










1.4 - 


2.1 


4.2 • 10~ 3 


8 


15 










2.1 - 


2.8 


1.5 • 10~ 3 


13 


18 






5.0 


10.0 


0.0 - 


0.7 


5.7 • 10~ 3 


7 


14 










0.7 - 


1.4 


4.9 • 10~ 3 


8 


17 










1.4 - 


2.1 


2.4 • 10~ 3 


12 


25 


10.0 


- 30.0 


1.7 


3.0 


0.0 - 


0.7 


1.7 • 10 ^ 


5 


13 










0.7 - 


1.4 


1.3 • 10~ 2 


6 


11 










1.4 - 


2.1 


9.8 • 10~ 3 


7 


10 










2.1 - 


2.8 


4.2- 10~ 3 


11 


14 










2.8 - 


3.5 


1.5 • 10~ 3 


20 


30 






3.0 


5.0 


0.0 - 


0.7 


1.5 • 10~ 2 


4 


12 










0.7 - 


1.4 


1.2- lO- 2 


4 


11 










1.4 - 


2.1 


8.8 • 10~ 3 


5 


12 










2.1 - 


2.8 


5.8 • 10~ 3 


7 


11 










2.8 - 


3.5 


1.2 • 10~ 3 


13 


13 






5.0 


1.0 


0.0 - 


0.7 


1.4- 10~ 2 


3 


12 










0.7 - 


1.4 


1.1 ■ 10~ 2 


3 


12 










1.4 - 


2.1 


7.5 • 10~ 3 


5 


15 










2.1 - 


2.8 


3.7 • 10~ 3 


7 


14 










2.8 - 


3.5 


8.1 • 10~ 4 


16 


30 






1.0 


33.0 


0.0 - 


0.7 


1.2- 10~ 2 


4 


13 










0.7 - 


1.4 


9.6 • 10~ 3 


4 


15 










1.4 - 


2.1 


5.2- 10~ 3 


7 


18 










2.1 - 


2.8 


1.2 ■ 10~ 3 


13 


28 


30.0 


- 100.0 


5.0 


10.0 


0.0 - 


0.7 


2.6 • 10- 7 


6 


11 










0.7 - 


1.4 


2.6 • 10~ 2 


6 


14 










1.4 - 


2.1 


1.6 • 10 2 


7 


12 










2.1 - 


2.8 


7.9 • 10~ 3 


11 


13 










2.8 - 


3.5 


3.8 • 10~ 3 


21 


24 






10.0 


- 100.0 


0.0 - 


0.7 


2.7 • 10~ 2 


3 


13 










0.7 - 


1.4 


2.2- 10~ 2 


3 


11 










1.4 - 


2.1 


1.4- 10~ 2 


4 


12 










2.1 - 


2.8 


4.9 • 10~ 3 


6 


14 










2.8 - 


3.5 


7.3 • 10~ 4 


15 


22 



Table 8: Dijet rate R 2 averaged over the regions indicated in x, Q 2 and \Arj*\ for A = 2 GeV. 
The measurement is restricted to values of the inelasticity variable y between 0.1 < y < 0.7 and 
to values of the polar angle of the scattered electron between 156° < 6 < 175°. 
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Q 2 

[GeV 2 ] 


x/10- 4 


S 


Ostat &syst. 
[%] [%] 


5.0 - 10.0 


1.0 - 1.7 
1.7 - 3.0 
3.0 - 5.0 
5.0 - 10.0 


0.086 
0.053 
0.049 
0.038 


8 9 

9 9 
13 9 
21 11 


10.0 - 15.0 


1.7 - 3.0 
3.0 - 5.0 
5.0 - 10.0 
10.0 - 18.0 


0.056 
0.048 
0.039 
0.022 


11 9 
13 12 
15 10 

33 24 


15.0 - 20.0 


3.0 - 5.0 

c f\ i a A 

5.0 - 10.0 
10.0 - 22.0 


0.066 

r\ ACA 

0.050 
0.039 


13 7 
15 16 
23 24 


20.0 - 30.0 


3.0 - 5.0 
5.0 - 10.0 
10.0 - 33.0 


0.086 
0.051 
0.043 


16 14 
12 12 
15 13 


30.0 - 50.0 


5.0 - 10.0 
10.0 - 55.0 


0.058 
0.038 


13 10 
11 16 


50.0 - 100.0 


10.0 - 25.0 
25.0 - 100.0 


0.040 
0.038 


14 16 
13 11 



Table 9: Measured ratio S for jets with an azimuthal separation of A0* < 120° for A = 2 GeV 
as shown in Figures|9]to[TT] The measurement is restricted to values of the inelasticity variable 
y between 0.1 < y < 0.7 and to values of the polar angle of the scattered electron between 
156° < 6 < 175°. 
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